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ABSTRACT 

As part of a large study to investigate the nature of the longer wavelength continuum emission 
of radio-loud AGN, we present new mid to far-infrared (MFIR) and high frequency radio obser- 
vations for a complete sample of 2Jy powerful, southern radio galaxies at intermediate redshifts 
(0.05 < z < 0.7). Utilizing the sensitivity of the Spitzer Space Telescope, we have made deep 
MIPS observations at the wavelengths of 24, 70 and 160/im, detecting 100% of our sample at 
24/im and 90% at 70/im. This high detection rate at MFIR wavelengths is unparallelled in sam- 
ples of intermediate redshift radio galaxies. Complementing these results, we also present new 
high frequency observations (15 to 24Ghz) from the Australia Telescope Compact Array (ATCA) 
and the Very Large Array ( VLA), which are used to investigate the potential contamination of the 
MFIR continuum by non-thermal synchrotron emission. With the latter observations we detect 
the compact cores in 59% of our complete sample, and deduce that non-thermal contamination 
of the MFIR continuum is significant in a maximum of 30% of our total sample. MFIR fluxes, 
radio fluxes and spectral energy distributions for the complete sample are presented here, while 
in a second paper we will analyse these data and discuss the implications for our understanding 
of the heating mechanism for the warm/cool dust, star formation in the host galaxies, and the 
unified schemes for powerful radio sources. 

Subject headings: galaxies:active-infrared:galaxies-radio:galaxies 



1. Introduction 

The mid to far-infrared (MFIR) is one of the 
most profitable wavelength ranges for investigat- 
ing AGN and associated phenomena. This is pri- 
marily because AGN emission suffers less from ef- 
fects of obscuration at these wavelengths. More- 
over, the thermal infrared samples the emission 
that is absorbed and then re-radiated by dust, in- 
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directly sampling the radiative power of the active 
core and any starburst component that may be 
present. Therefore this wavelength range is cru- 
cial both for testin g the unified schemes for AGN 
(e.g. lBarthe]|[l989h . and investigating the link be- 
tween AGN and starburst activ ity as part of the 
evolution of the host galaxies (iRowan-Robinson 
19951 lHaas et al.lll998l [Archibald et al.ll200ll ). 



In terms of testing unified schemes, radio-loud 
samples of AGN have the advantage that their ex- 
tended radio emission is isotropic, allowing sam- 
ples to be selected free from the orientation bias 
that affects samples selected at X-ray, optical and 
near-IR wavelengths. However, in order to fully 
realize the potential of MFIR observations for 
studies of radio galaxies it is important to obtain 
observations of carefully selected samples that at- 
tain a high level of completeness in terms of de- 
tections. Unfortunately, observing in the MFIR 
is technically more challenging than at near- 
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infrared wavelengths. Therefore, although the 
IRAS observatory made many significant contri- 
butions to our understanding of the MFIR infrared 
prop e rties of radio- loud AGN dNeugebauer et alj 
19861 iGolombek et alj 1 19881 iKnapp et al l | l99oT 
Impev fc Gregqrini 19931 Irleckman et al.l Il992l 
1994l lHes et al.lll99^ . it had a relatively low sen- 
sitivity, resulting in a typical MFIR detection rate 
of less than 30%, even for radio galaxies in the 
local universe (z < 0.3). 

Subsequent attempts were made to observe 
samp les of radio galaxies w i th the ISO observa- 
tory dHaas et al J 1 19981 12004 IPolletta et alJ fe oOO 



van Bemmel et al.l2000LlMeisenheimer et al. 2001 ) 
However, meager, sometimes heterogeneous sam- 
ples, and the modest sensitivity of the observa- 
tory, only allowed for a small improvement on 
the achievements of IRAS. Studies of samples of 
radio-loud AGN with ISO have typical detection 
rates of no more than 50% in the far-IR, impeding 
attempts to understand the MFIR emission. 

The launch of Spitzer in 2003, brought an op- 
portunity to make observations with orders of 
magnitude improved sensitivity at MFIR wave- 
lengths. However, published Spitzer studies of ra- 
dio loud AG N are either bas ed on heterogeneous 
samples (e.g. IShi et al" 2005 ). or on samples that 
are relatively faint and distant . resulting in a low - 
detections rates at 7 0/xm (e.g. IClearv et al.ll2007 , 
Seymour et"aTll2007l ). 

A further problem is that many previous MFIR 
studies of radio-loud AGN samples are selected 
from the 3C radio catalogue, which are ham- 
pered by the lack of published high quality op- 
tical spectroscopic observations for many of the 
objects. In contrast, the southern 2Jy sample 



(|Tadhunter et all Il993h 



is 



unique in the sense 
that deep spe ctra have been published for the 



Wills et al 



, , Pi 

whole sample (ITadhunter et al.lll993L 1 19981 liool 



2002, 



Holt et al. 200 



Measured 

emission line luminosities for this sample pro- 
vide in formation about the in trinsic power of the 
AGN (jTadhunter et al.l Il998l ). Moreover, 



care- 



ful modelling of the continuum spectra provides 
key information about the stellar populations in 
the host galaxies and, in partic ular, in the pres- 



ence of young stellar po pulations (jTadhunter et al 



2002; I Wills et al.l 12004 ). The 2Jy sample has also 



been extensively obs e rved at rad io wavelengths 
(|Morganti et al.lll993l Il997l 1 19991 ). allowing links 



to be made between radio and optical properties 
and independent estimates of the orientation of 
the sources to the line of sight. In view of its com- 
pleteness and the availability of deep spectroscopic 
and radio data, which makes it well suited to in- 
vestigating the nature of the MFIR emission and 
testing the unified schemes, we have undertaken 
a program of deep imaging with Spitzer /MIPS of 
the 2Jy sample. 

Given that we are interested in the thermal 
emission and heating mechanism of the warm/cool 
dust sampled through the MFIR fluxes, possi- 
ble contamination from non-thermal synchrotron 
sources could potentially bias our conclusion. 
Therefore, to complement the MFIR observations, 
we have also undertaken high frequency observa- 
tions at radio wavelengths (15 to 22 GHz) using 
ATCA and the VLA, with the aim of investigating 
the potential degree of non-thermal contamina- 
tion. 

A preliminary analysis of the results from thi s 
program was presented in ITadhunter et al~ ( 2007 ). 
Here we present the MFIR and radio core data, 
and discuss the potential non-thermal contamina- 
tion of the MFIR continuum. A second paper com- 
prises an in-depth analysis of these data (Dicken 
et al. 2008, in preparation). 

2. Sample selection 

The sample for this study comprises a com- 
plete sub-sample of all 46 steep-spectrum pow- 
erful radio galaxies and quasars with redshifts 
0.05 < z < 0.7, and fl ux densities S^.jqnr. > 
2 Jy from the sample of ITadhunter et alJ (|l993l ) . 
We define the the steep spectrum selection as 



'2.7 



> 0.5 {F v oc v- 



which excludes 16 ob- 



jects within the redshift range. We also include 
PKS0347+05, which has since proved to fulfil the 



same selection criteria (|di Serego-Alighieri et al 
19941 ) . The lower redshift limit has been set to 



ensure that these galaxies are genuinely powerful 
sources, and the steep spectrum selection for the 
quasars rules out objects dominated by emission 
from the beamed relativistic jet and core compo- 
nents. Overall, the full sample of 46 objects in- 
cluded a mixture of broad line radio galaxies/radio 
loud quasars (BLRG/Q: 35%), narrow line radio 
galaxies (NLRG: 43%) and weak line radio galax- 
ies (WLRG: 22%). In terms of radio morphology 
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classification, our complete sample includes 72% 
FRII sources, 13% FRI sources and 15% compact 
steep spectrum/gigahertz peak spectrum objects. 

For comparison purposes we also observed or 
collected data from the Spitzer/IRAS archive, 
for the flat spectrum, core dominated objects 
3C273, PKS0521-36 and PKS1549-79. Note 
that PKS1549— 79 is a particularly interesting 
source, because it a rare example of a powerful, flat 
spectrum radio source that is classified as a galaxy 
at optical wavelength s; the nature of this object is 
discussed in detail in lHolt et al. ( 20061 ). However, 



these objects are not part of what we refer to as 
our complete sample. Our complete sample of 46 
objects comprises the 49 objects with Spitzer or 
IRAS observations presented below, minus the 3 
comparison flat spectrum objects. 

3. MFIR observations and reduction 

3.1. Observations 

We have obtained new observations for 43 
objects using the Spitzer Space Telescope and, 
for 37 of these, we have obtained the very first 
MFIR detections. In addition, we include re- 
sults from data obtained for 5 other targets from 
the Spitzer archive Reserved Observations Cata- 
logue (ROC): PKS0915-11 (3C218), PKS1226+02 
(3C273), PKS1559+02 (3C327), PKS1648+05 
(3C348) and PKS2221-02 (3C445). The obser- 
vations were carried out between August 2005 
and Jan 2007, and the ROC targets were observed 
between March and November 2004. 

All targets were observed with the Multi- 
band Imag i ng P hotometer for Spitzer (MIPS: 
Rieke et al.l 120041) at wavelengths of 24, 70 and 
160/im, apart from 5 targets that were not ob- 
served at 160/im due to scheduling reasons. De- 
tails of the observations are shown in Table [T] 

3.2. Reduction 

The MIPS instrument on board Spitzer takes 
short exposure images in a dither pattern. These 
are then combined to make a final 'mosaiked' 
image. Rather than use the pipeline processed 
data, we have re-reduced data for all the 48 ob- 
jects in the sample starting from the Basic Cal- 
ibrated Data (BCD) files to produce final mo- 
saiked images for analysis. The reduction was car- 



ried out at Rochester Institute of Technology and 
The University of Sheffi eld using the MOPEX soft- 
ware reduction package ( M akovoz et al . 2006) pro- 
vided by the Spitzer Science Center (SSC). Due to 
the different nature of the detectors and sensitiv- 
ities at the 3 wavelengths observed, the reduction 
method we have selected is different for each. Thus 
we discuss the 3 bands separately below. 

3.3. 24/im Reduction 

As the shortest wavelength band, the 24/im 
channel of MIPS is also the most sensitive, and 
has the highest spatial resolution (6" FWHM). In- 
tegration times for our targets were in the range 
48.2 to 180.4 seconds, where over 80% have the 
longer exposure time (see Table [T]). 

Minimal processing to the 24/im SSC pipeline 
product is required for most general science ob- 
jectives. However, we re-processed the data by 
mosaiking the BCD files using the MOPEX soft- 
ware and adding an additional flatfield and over- 
lap correction. The flatfield script is part of the 
MOPEX software package. It computes a flat field 
from the median of all the dithered frames, nor- 
malizes this to an average of 1 and divides into all 
the input BCD files. In contrast, the overlap cor- 
rection interpolates the input images onto a com- 
mon grid and then the cumulative pixel-by-pixel 
difference between the overlapping areas is mini- 
mized. The final mosaic pixel size was set to 2.45 
arsecs. This was deliberately chosen to match the 
pipeline post-BCD products for purposes of com- 
parison. Examples of the images can be seen in 
Figure Q] 

The 24/im fluxes were extracted using aper- 
ture photometry in the GAIA package. Aper- 
ture corrections were derived from an empirically 
determined average curve of growth of flux vs. 
aperture diameter derived from observations of 
bright sources in our sample with S24 > 50mJy 
(PKS0521-36, PKS1226+02, PKS1949+02, PKS2221- 
02, PKS2314+03). These aperture corrections are 
broadly consistent with those published by the 
SSC. For most objects the aperture was set to a 
standard radius of 15" corresponding to an aper- 
ture correction of 1.08. However, for a few objects 
an aperture of half this size (7.5") was used to 
avoid contamination from the flux of objects close 
to the source. In the latter cases the applied aper- 
ture correction was 1.52. 
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Table 1 

Spitzer Sample and Observational Data 



PKS Name 


Other 


RA(J2000) 


Dcc(J2000) 


Opt. Class 


Rad. Class 


z 


Obs. Date ti 


„ t 24(s) 


tint70(s) 


tintl60(l 


DO'?'} — 9fi 




00 


25 


49.18 


-26 


02 


12 


.8 


NLRG 


CSS 


322 


Vlcr 05 


180 4 


543 


3 


167 8 


0034 — 01 


3C015 


00 


37 


04.14 


-01 


09 


08. 


2 


WLRG 


FRII 


0.073 


Jan-06 


180.4 


543 


.3 


167.8 


0035 — 02 


3C17 


00 


38 


20.51 


-02 


07 


40. 


1 


BLRG 


(FRII) 


0.220 


Jul-06 


180.4 


543 


.3 


167.8 




3C18 


00 


40 


50.53 


+ 10 


03 


26. 


.8 


blrg 


FRII 


0. 188 


Jan-06 


180.4 


543 


.3 


167.8 


nnsQ— 44 




00 


42 


09.03 


-44 


14 


01. 


3 


NLRG 


FRII 


0.346 


Dcc-05 


180.4 


543 


3 


167.8 


0043 — 42 




00 


46 


17.75 


-42 


07 


51. 


4 


WLRG 


FRII 


0.116 


Dec-05 


180.4 


543 


.3 


167.8 


0105 — 16 


3C32 


01 


08 


16.90 


-16 


04 


20. 


6 


NLRG 


FRII 


0.400 


Jul-06 


180.4 


543 


.3 


167.8 


0117— 15 


3C38 


01 


20 


27.16 


-15 


20 


17. 


8 


NLRG 


FRII 


0.565 


Jul-06 


180.4 


543 


.3 


167.8 


0213 — 13 


3C62 


02 


15 


37.50 


-12 


59 


30. 


.5 


NLRG 


FRII 


0.147 


Jan-07 


180.4 


545 


.3 


167.8 


0235 — 1 9 


OD-159 


02 


37 


43.45 


-19 


32 


33. 


3 


BLRG 


FRII 


0.620 


Feb- 06 


180.4 


545 


.3 


167.8 


0252—71 




02 


52 


46.26 


-71 


04 


35. 


9 


NLRG 


CSS 


0.566 


Aug-05 


180.4 


545. 


.3 


167.8 


03474-05 




03 


40 


46.50 


+05 


51 


42. 


3 


BLRG 


FRII 


0.339 


Feb- 06 


180.4 


545 


.3 


167.8 


0349—27 




03 


51 


35.81 


-27 


44 


33 


.8 


NLRG 


FRII 


0.066 


Aug-05 


180.4 


545 


3 


167.8 


04044-03 


3C105 


04 


07 


16.49 


+03 


42 


25. 


.8 


NLRG 


FRII 


0.089 


Sen-05 


180.4 


545 


.3 




0409 — 75 




04 


08 


48.49 


-75 


07 


19. 


3 


NLRG 


FRII 


0.693 


Aug-05 


180.4 


545 


3 


167.8 


0442 — 28 




04 


44 


37.67 


-28 


09 


54. 


6 


NLRG 


FRII 


0.147 


Fcb-06 


180.4 


545 


.3 


167.8 


0521 —36 

Y.) \J _L -.i V.I 




05 


22 


58.94 


-36 


27 


31. 


2 


BLRG 


C/J 


0.055 


Dcc-05 


92.3 


230 


.7 


167.8 


0620—52 

UUi W Oil 




06 


21 


43.29 


-52 


41 


33 


3 




FRI 


0.051 


Dcc-05 


180.4 


545 


3 


167.8 


0625 — 35 


OH-342 


06 


27 


06.65 


-35 


29 


16. 


3 


WLRG 


FRI 


0.055 


Dcc-05 


180.4 


545 


.3 


167.8 


0625—53 




06 


26 


20.44 


-53 


41 


35. 


2 


WLRG 


FRII 


0.054 


Dcc-05 


180.4 


545 


3 


167.8 


0806 — 10 


3C195 


08 


08 


53.64 


-10 


27 


40. 


2 


NLRG 


FRII 


0.110 


Dec-05 


180.4 


545 


3 


167.8 


0859 — 25 




09 


01 


47.50 


-25 


55 


19. 





NLRG 


FRII 


0.305 


Dcc-05 


180.4 


545 


3 


167.8 


0915 — 11 


Hydra A 09 


18 


05.67 


-12 


05 


44. 





WLRG 


FRI 


0.054 


May-04 


165.7 


125 


.8 


83.9 


0945+07 


3C227 


09 


47 


45.15 


+07 


25 


20. 


4 


BLRG 


FRII 


0.086 


Dec-05 


180.4 


545 


.3 


167.8 


1136-13 




11 


39 


10.75 


-13 


50 


43 


.1 


Q 


FRII 


0.554 


Jun-06 


180.4 


545 


3 


167.8 


1151-34 




11 


54 


21.79 


-35 


05 


29 


.1 


Q 


CSS 


0.258 


Feb-06 


180.4 


543 


.3 


167.8 


1226+02 


3C273 


12 


20 


06.70 


+02 


03 


08. 


6 


Q 


C/J 


0.158 


Jun-04 


48.2 


37 


.7 


41.9 


1306-09 




13 


08 


39.17 


-09 


50 


32 





NLRG 


CSS 


0.464 


Feb-06 


180.4 


545 


3 


167.8 


1355-41 




13 


59 


00.23 


-41 


52 


54 


.1 


Q 


FRII 


0.313 


Feb-06 


180.4 


545 


.3 


167.8 


1547-79 




15 


55 


21.66 


-79 


40 


36. 


3 


BLRG 


FRII 


0.483 


Aug-05 


180.4 


545 


3 


167.8 


1559+02 


3C327 


16 


02 


27.38 


+01 


57 


55. 


.7 


NLRG 


FRII 


0.104 


Mar-04 


165.7 


125 


.8 


83.9 


1602+01 


3C327.1 


16 


04 


45.35 


+01 


17 


50. 


.8 


BLRG 


FRII 


0.462 


Aug-05 


180.4 


545 


.3 


167.8 


1648+05 


Here A 


16 


51 


08.16 


+04 


59 


33. 


8 


WLRG 


FRI 


0.154 


Mar-04 


165.7 


125 


.8 


83.9 


1733-56 




17 


37 


35.80 


-56 


34 


03 


4 


BLRG 


FRII 


0.098 


Aug-05 


180.4 


545 


.3 


167.8 


1814-63 




18 


19 


34.96 


-63 


45 


48. 


.1 


NLRG 


CSS 


0.063 


Aug-05 


180.4 


545 


3 


167.8 


1839-48 




18 


43 


14.64 


-48 


36 


23 


3 


WLRG 


FRI 


0.112 


Oct-06 


180.4 


545 


3 


none 


1932-46 




19 


35 


56.56 


-46 


20 


40 


.7 


BLRG 


FRII 


0.231 


May-06 


180.4 


545 


3 


167.8 


1934-63 




19 


39 


24.99 


-63 


42 


45. 


6 


NLRG 


GPS 


0.183 


May-06 


180.4 


545 


3 


167.8 


1938-15 




19 


41 


15.15 


-15 


24 


30. 


9 


BLRG 


FRII 


0.452 


Nov-06 


180.4 


545 


.3 


none 


1949+02 


3C403 


19 


52 


15.77 


+02 


30 


23. 


1 


NLRG 


FRII 


0.059 


Nov-05 


92.3 


230 


.7 


167.8 


1954-55 




19 


58 


16.06 


-55 


09 


37 


.5 


WLRG 


FRI 


0.060 


May-06 


180.4 


545 


3 


167.8 


2135-14 




21 


37 


45.18 


-14 


32 


55. 


.5 


Q 


FRII 


0.200 


Nov-05 


180.4 


545 


.3 


167.8 


2135-20 


OX-258 


21 


37 


50.00 


-20 


42 


31. 


.7 


BLRG 


CSS 


0.635 


Nov-05 


180.4 


545 


.3 


167.8 


2211-17 


3C444 


22 


14 


25.76 


-17 


01 


36. 


2 


WLRG 


FRII 


0.153 


Dec-06 


180.4 


545 


.3 


none 


2221-02 


3C445 


22 


23 


49.57 


-02 


06 


12. 


4 


BLRG 


FRII 


0.057 


Nov-04 


92.3 


69 


2 


none 


2250-41 




22 


53 


03.16 


-40 


57 


46. 


2 


NLRG 


FRII 


0.310 


Nov-05 


180.4 


545 


3 


167.8 


2314+03 


3C459 


23 


16 


35.21 


+04 


05 


18. 


2 


NLRG 


FRII 


0.220 


Nov-05 


92.3 


230 


.7 


167.8 


2356-61 




23 


59 


04.50 


-60 


54 


59 


.1 


NLRG 


FRII 


0.096 


Dec-05 


180.4 


545 


3 


167.8 



Note. — Tab lc^ Column 5: broad line radio galaxy (BLRG), narrow line radio galaxy (NLRG), quasar (Q), weak line radio 
galaxy (WLRG, also known as low excitation galaxies). Column 6: Fan ar off- Riley class 1 (FRI), Fan aroff- Riley class 2 (FRII), 
compact steep spectrum type (CSS), core/jet (C/J), gigahertz peaked spectrum (GPS). Uncertain classification in brackets. 
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Table 2 
Spitzer Data Results 



PKS Name 


Other name 


z 


S24 (mjy) 


cr 


S70 (mjy) 


cr 


Sieo (mjy) 


a 






322 


2.4 


3 


26 3 


3 1 


83 6 


16 7 


nnid— m 


3C015 


073 


7.5 


2 


17 9 


2 1 


< 43 1 


14.4 




3C17 


220 


12 2 


1 


23 6 


4 8 


97 2 


9 1 


U W O O \J V 


3C18 


188 


25 9 


3 


32 2 


4.8 


< 37 8 


12 6 


uuoa — '-i L ± 




346 


33 


4 


68 7 


4 6 


84 6 


13 5 


0043—42 




0.116 


11.1 


0.2 


9.9 


3.0 


< 24.9 


8.3 


ninr;_ifi 


3C32 


400 


9.7 


2 


< 11 8 


3 9 


< 22 5 


7.5 


0117—15 


3C38 


0.565 


6.1 


0.2 


20.2 


1.8 


< 46.4 


15.5 


(191 3_1 ^ 


3C62 


0. 147 


40 2 


1 


37 1 


2 9 


< 51 1 


17.0 


f)9'i c i— 1 q 


OD-159 


620 


11 1 


2 


14 3 


2 5 


< 23 6 


7.9 


n 1 ?^? — 71 




5(56 


2 9 


1 


< 9 1 


3 


< 37.4 


12.5 






0.339 


3.5 


0.2 


30.8 


4.3 


< 52.2 


17.4 






066 


8.8 


3 


41 9 


1 3 


< 31 3 


10.4 




3C105 


089 


30 8 


1 


70 9 


2 










693 


1 5 


3 


11 2 


1 


< 36 4 


12 1 


0442 — 28 




0. 147 


22.0 


0.3 


31.0 


4.5 


< 25.9 


8.6 


C\ r \9'\ — ^fi 
yjo^i J. — ou 




Q55 


204 1 


3 


526 4 


6 


950 1 


13 3 






051 


4.5 


1 


47.3 


1.4 


< 24 3 


8.1 




OH-342 


055 


24.7 


3 


44.8 


1 3 


131 


11 2 


0fi95 — 53 




0.054 


1.7 


0.2 


< 10.8 


3 G 


< 52.6 


17.5 


UoUD — 1 u 


3C195 


110 


258 3 


4 


489 7 


1 3 


308 7 


7 1 


085Q — 95 




0.305 


9.3 


0.4 


8.4 


2.8 


< 52.3 


17.4 


0915 — 11 


Hydra A 


0.054 


8.9 


0.2 


115.3 


4.8 


164.4 


7.3 


0945+07 


3C227 


0.086 


47.7 


0.3 


19.4 


3.5 


< 44.1 


14.7 


1136-13 




0.554 


13.8 


0.2 


23.9 


2.8 


< 23.2 


7.7 


1151-34 




0.258 


16.4 


0.3 


51.5 


3.0 


< 46.3 


15.4 


1226+02 


3C273 


0.158 


499.1 


0.8 


414.9 


3.9 


402.1 


11.9 


1306-09 




0.464 


4.6 


0.2 


21.7 


2.1 


< 29.2 


9.7 


1355-41 




0.313 


53.1 


0.3 


66.3 


2.0 


< 55.9 


18.6 


1547-79 




0.483 


7.9 


0.1 


19.2 


1.7 


< 50.2 


16.7 


1559+02 


3C327 


0.104 


242.0 


0.3 


470.0 


3.7 


262.8 


21.4 


1602+01 


3C327.1 


0.462 


7.7 


0.3 


12.3 


2.5 


< 32.4 


10.8 


1648+05 


Here A 


0.154 


2.0 


0.2 


< 18.8 


6.3 


134.3 


15.8 


1733-56 




0.098 


29.2 


0.3 


151.0 


3.3 


318.2 


10.9 


1814-63 




0.063 


60.6 


0.4 


142.1 


2.1 


202.1 


13.1 


1839-48 




0.112 


3.1 


0.3 


10.9 


2.5 






1932-46 




0.231 


2.5 


0.1 


17.6 


1.6 


< 36.5 


12.2 


1934-63 




0.183 


17.4 


0.1 


19.9 


2.2 


< 33.2 


11.1 


1938-15 




0.452 


6.8 


0.3 


19.7 


4.2 






1949+02 


3C403 


0.059 


193.0 


0.2 


348.4 


3.7 


251.3 


7.1 


1954-55 




0.060 


2.7 


0.2 


8.8 


2.9 


< 34.1 


11.4 


2135-14 




0.200 


104.9 


0.2 


113.7 


4.8 


118.6 


20.5 


2135-20 


OX-258 


0.635 


4.3 


0.3 


37.0 


4.4 


127.1 


10.3 


2211-17 


3C444 


0.153 


0.5 


0.1 


< 9.4 


3.1 






2221-02 


3C445 


0.057 


232.1 


0.3 


186.4 


5.2 






2250-41 




0.310 


11.6 


0.1 


22.0 


2.5 


< 27.6 


9.2 


2314+03 


3C459 


0.220 


49.9 


0.3 


511.5 


3.9 


490.4 


19.1 


2356-61 




0.096 


41.0 


0.2 


74.7 


2.2 


< 35.4 


11.8 



Note. — Tablj2] All images reduced were from the BCD files, Spitzer Science Center pipeline 
version 14.4.0, Fluxes presented for 2135-20 are the mean fluxes of two identical observations see 

ED 
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Fig. 1. — Example of a MIPS mosaiked images of PKS1814— 63. Left: 24/iin image with an example aperture. 
Top Right: 70/im image with an example aperture. Bottom Right: 160/im image. The white circles indicate 
the apertures used for the object flux measurements, while the yellow circles delimit the annuli used to 
estimate the background sky flux. The bar in the top right hand corner of each image is of length 1 arcmin. 



The standard deviation of 6 aperture flux mea- 
surements of background sky patches, obtained us- 
ing an identical aperture to the object flux mea- 
surements, was used to derive the la uncertain- 
ties presented in Table [2] These la values repre- 
sent the fluctuations in the sky background due 
to photon counting noise, as well as mosaiking 
and flat fielding artifacts. Given that most of the 
sources are faint relative to the background, these 
la measurements give a realistic estimate of the 
la uncertainties in the flux measurement for the 
fainter sources in the sample. In addition there 
i s an estimated ±4% flu x calibration uncertainty 
(jEngelbracht et al.ll2007l ). which is likely to domi- 
nate for the brightest sources in the sample. Over- 
all, at 24/im we detect 100% of our sample at the 
> 5a level. 

3.4. 70/j,m Reduction 

The 70/im waveband is less sensitive than the 
24/im band and uses an entirely different Ge:Ga 
detector technology, which is not as stable as the 
Si:As detectors used at 24/im. The spatial reso- 
lution of Spitzer/MIPS at 70/im is 18"(FWHM). 
Integration times for our targets were in the range 



37.7 to 545.3 seconds, where over 80% have the 
longer exposure time (see Table [T]). The 70/tm 
data were also reduced with the MOPEX pack- 
age. In addition to the basic mosaiking process we 
used a column filtering process from the SSC con- 
tributed software pages (IDL program BCD col- 
umn filter. pro), which calculates the median value 
for each column and subtracts that from each col- 
umn for each BCD. This process significantly im- 
proves artifacts in the images with minimal loss of 
flux from the source. The final mosaic pixel size 
is 3.9", again chosen to match that of the pipeline 
mosaic products. 

The fluxes at 70/im were also measured us- 
ing aperture photometry, with an aperture ra- 
dius of 24" for the majority of sources. Aper- 
ture corrections and la uncertainties were derived 
in the same manner as for the 24/im imagetQ. 
Again these aperture corrections are broadly con- 



X N.B. In the case of 70/im and 160/im, in addition to the 
photon counting and instrumental induced fluctuations, 
there will be an additional contribution to the fluctuations 
in the background from sources close to the detection limit. 
This will be reflected in the lcr estimates tabulated in Table 

m 
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sistent with those published by the SSC. In this 
case the aperture correction factor is 1.33. A 
smaller aperture of half the size (22") was used 
for fainter sources to obtain better S/N for de- 
tections, with a correction factor of 2.04. There is 
also an additional ±1 % calibration error at 70/im 
([Gordon et al.1 120071 ). For undetected sources up- 
per limits were derived using the standard devia- 
tion of background measurements to obtain a. For 
these sources the values presented in Table [5] arc 
the 3<T upper limits, of which there are 5. 

Overall, we detect 90% of our sample at 70/im 
at the > 3tr level. 

3.5. 160/im Reduction 

The longest wavelength MIPS band (160/im) 
has the lowest effective sensitivity and utilizes 
the same detector technology as 70/im. The 
spatial resolution of Spitzer/ MIPS at 160/im is 
40"(FWHM). Integration times for our targets at 
160/im were in the range 41.9 to 167.8 seconds, 
where over 90% have the longer exposure time 
(see Table CD). 

In this case we mosaiked the images with the 
MOPEX software (final image pixel size 7.9") with 
no additional processing, as none of the available 
tools appeared to provide any significant improve- 
ment to the images. 

Aperture photometry proved impossible to use 
for extracting the 160/im flux, due to the dif- 
ficulty in estimating the background accurately. 
This is because the PSF at 160/im is on the scale 
of the image, and there is also the strong pos- 
sibility of contamination by other sources in the 
field. Therefore, we used the APEX program in 
the MOPEX GUI software to extract the fluxes 
of detected sources using PSF fitting. The stan- 
dard PSF available with the software was used for 
this task. Overall we detect 33% of our sample at 
160/mi at the > 3a level. 

The lcr errors were derived in an identical way 
to the other 2 wavebands, and the upper limits 
were derived in an identical way to 70/im . There 
is also an additional ±20% flux calibration uncer- 
tainty as stated in the MIPS data handbook (avail- 
able from the SSC website). Upper limits were de- 
rived using the standard deviation of background 
measurements to obtain a The values presented 
in Tableware the 3<r upper limits, of which there 



are 27. Note that a few objects undetected by the 
APEX software appear to be detected in a visual 
inspection of the image. Also, overlapping PSF 
from nearby objects has prevented the APEX de- 
tection of a point source in at least one image. 
Therefore, we believe that the APEX software un- 
derestimates the true detection rate of sources in 
our sample; in this sense the detection rate is con- 
servative. This conclusion is consistent with the 
fact that the typical 3cr detection limit for the 
160/tm data is 55mJy, which is somewhat lower 
than the lowest measured flux of 84.6mJy. 

3.6. Notes on the MFIR observations 

Due to an error in the acquisition file for the 
data, PKS2135-20 was observed twice. This gave 
us an opportunity to investigate the uncertainty 
derived from two identical observations of the 
same source. We found that the fluxes measured 
for this source in the two data sets were consistent 
within 10% in all MIPS bands. Fluxes presented 
in Table 2 are the mean fluxes of the two observa- 
tions. 

It has also been possible to compare some of 
our flux values obtained using Spitzer with those 
based on IRAS observations. Taking into account 
BLRG/quasar objects that are potentially vari- 
able, we directly compared the IRAS 25//m and 
Spitzer 24/tm flux measurement. However in order 
to make the comparison for the FIR fluxes it was 
necessary to extrapolate the IRAS 60/im fluxes to 
the Spitzer 70/im band using the IRAS 25 to 60/im 
spectral index. We find that the Spitzer and IRAS 
fluxes agree to within better than 20% at 24/im 
and within better than 30% at 70/im. It is also 
worth noting that our Spitzer measured flux val- 
ues for sources i n commo n agre e with published 
Spitzer fluxes in IShi et al.1 ( 2005 ) to within a few 
percent. 

4. Radio Observations and Reduction 
4.1. Sample 

The radio sample is identical to that discussed 
in section[5J The objects were divided between the 
ATCA and the VLA observatories according to the 
declination of the source (cutoff 8 » —25°). Four 
objects - PKS0023-26, PKS0349-27, PKS0442-28, 
PKS0859-25- were observed with both ATCA and 
the VLA. 
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Fig. 2. — Example images from our high frequency radio observations. Top: PKS1355— 41, an FRII quasar. 
Bottom: PKS1932— 46, an FRII broad line radio galaxy. Contour level for the plots are as follows: 2.4*-l, 1, 
2, 4, 6, 8, 12, 16, 32, 64, 128 mJy/beam. The core and the lobes/hot spots are clearly seen in both objects. 
In PKS 1932-46, residuals are seen on large scales due to extended, low surface brightness structures that 
cannot be properly imaged due to the poor uv coverage. 



Table 3 
ATCA Radio Core Data 



PKS Name 


Radio Class 


z 


S co™ ( mJ y) 


a 


Bcam(" ) 


Score ( m Jy) 


<7 


Bcam( ) 


Combined R 


nmo aa 


FRII 




<, 1.1 


n i 


n ^wn *zf 9 7"\ 
u.oxu.o^i ( ) 


n q 
<^ u.y 


U.O 


n /fvn oi o 

U.4XU.Z 1 -£i, o ) 


^n nn^o 


nri/i i ao 


FRII 


U . X X o 




U.o 


n Rvfl ^f OR /l ^ 




U.o 


U.OXU.Z^O-L.O J 


1.4 U.UUUO 


flOCO 71 




0.566 


278. 1 


1.0 


U.4XU.OI lo.^ J 


187.7 


1.0 


U.oxu. z ^-zu.oj 


^n nifin 


Uo4y — ^ / 


T7D TT 


0.066 


20.0 


0. 1 


U. /xu.o^o.o j 


21.7 


0.3 


U.OXU.^4. / ) 


u.uuyo 


0409-75 


FRII 


0.693 


< 10.0 


3.3 


0.4x0.3(24.8) 


< 9.0 


3.0 


0.3x0.2(21.6) 


<0.0023 


0442-28 


FRII 


0.147 


19.1 


0.4 


0.9x0.3(19.2) 


20.4 


0.4 


0.6x0.2(17.4) 


0.0090 


0521-36 


C/J 


0.055 


2798* 


24.4 


0.6x0.4(11.5) 


2654* 


14.8 


0.4x0.3(14.5) 


0.4115 


0620-52 


FRI 


0.051 


118.6 


0.6 


0.5x0.3(26.3) 


101.2 


0.7 


0.4x0.2(28.6) 


0.0964 


0625-35 


FRI 


0.055 


571.6 


3.4 


0.7x0.3(6.2) 


478.8 


3.0 


0.5x0.2(7.5) 


0.3293 


0625-53 


FRII 


0.054 


20.4 


0.2 


0.5x0.2(7.0) 


19.3 


0.2 


0.3x0.2(8.6) 


0.0130 


0859-25 


FRII 


0.305 




0.5 


1.5x0.3(2.2) 




0.5 


1.1x0.2(3.6) 


2.5 <0.0013 


1151-34 


CSS 


0.258 


759.9* 


6.5 


1.0x0.3(-23.4) 


598.2* 


6.3 


0.7x0.2(-22.0) 


<0.0072 


1355-41 


FRII 


0.313 


102.7 


0.8 


0.5x0.3(0.2) 


89.6 


1.0 


0.4x0.3(-0.8) 


0.0715 


1547-79 


FRII 


0.483 




0.2 


0.6x0.3(-18.4) 




0.2 


0.5x0.2(-14.8) 


1.4 0.0010 


1733-56 


FRII 


0.098 


293.6 


0.9 


0.4x0.3(-34.1) 


269.4 


1.9 


0.3x0.2(-32.5) 


0.0911 


1814-63 


CSS 


0.063 


1909* 


6.5 


0.4x0.3(-5.3) 


1626* 


7.3 


0.3x0.2(-5.0) 


0.0265 


1839-48 


FRI 


0.112 


108.2 


0.5 


0.5x0.3(-3.3) 


100.9 


0.6 


0.4x0.2(-3.3) 


0.0889 


1932-46 


FRII 


0.231 


16.4 


0.7 


0.5x0.4(-67.4) 


21.5 


0.5 


0.4x0.3(-65) 


0.0050 


1934-63 


GPS 


0.183 


1045* 


2.0 


0.6x0.3(68.9) 


755.1* 


1.1 


0.6x0.2(66.3) 




1954-55 


FRI 


0.060 


61.7 


0.2 


0.4x0.3(56.6) 


61.6 


0.3 


0.3x0.3(63.8) 


0.0393 


2250-41 


FRII 


0.310 


< 1.9 


0.4 


0.7x0.3(-20.2) 


< 1.1 


0.4 


0.5x0.3(-19.2) 


<0.0012 


2356-61 


FRII 


0.096 


59.7 


0.3 


0.4x0.3(76.6) 


63.8 


0.4 


0.3x0.2(89.4) 


0.0137 



Note. — Tablc[3] Flux density results from the ATCA observations. Definitions for column 2 see Tabltf^ (*) denotes a total flux 
measurement as the resolution of the observations is not high enough to resolve the core in these compact objects. Column 11 gives the 
R parameter defined as S CO re/{Stot — S cor e) where the total flux is taken at 5GHz assuming a fiat spectrum. 



4.2. ATCA Observations 

We observed 23 objects (with 8 < -25°) 
from our sample with ATCA. For one object 
(PKS0023-26) the quality of the ATCA data 
turned out to be quite low, we have therefore 
excluded it from our ATCA list and we will use 
only the VLA data. 

The sources were observed using the 6-km con- 
figuration - the longest available for this radio tele- 
scope. These observations were carried out be- 
tween the 19th and 23rd of July 2006 in perfect 
weather conditions. The total run comprised of a 
continuous lOOh of observations. 

Data were taken simultaneously at 17.9 GHz 
and 24.1 GHz. These high frequencies allow us to 
achieve a relatively high resolution. Each source 
was observed for a period ranging from 30 to 
60 min (depending on scheduling constraints) in 
scans of 10 min spread over a period of 12h in 
order to obtain a good enough uv-coverage. In 
each scan we swapped between the two frequen- 
cies. Furthermore, for each frequency, we actually 
took data simultaneously at 17.47 and 18.50 GHz 
in the first observation and at 23.62 and 24.64 GHz 



in the second, with a bandwidth of 128 MHz for 
each of these frequencies. These separate frequen- 
cies allowed us to improve (radially) the uv cover- 
age. The flux calibrator - PKS1934-64 flux level at 
12mm — 1.03Jy - was observed few times during 
the whole run. In addition, a nearby phase cali- 
brator was observed for 5 min before every scan 
on a source. A reference pointing observation was 
also carried out before every source. 

The ATCA data were cali brated us i ng the 



19951) . 



MIRIAD software package (jSault et al. 
The final images were obtained using multi- 
frequency synthesis option and adding together 
the two nearby frequencies (17.47 and 18.50 GHz; 
23.62 and 24.64 GHz) to increase the sensitiv- 
ity and improve the uv-coverage. Typically two 
cycles of self-calibration were needed to obtain 
the final images. In these images, we have used 
uniform weighting to achieve the highest possible 
resolution. 

Table [3] summarises the results of the ATCA ob- 
servations, including beam size and noise for each 
object and each frequency. We have detected radio 
cores in 14 of the 22 sources (63%; PKS 0023-26 is 
excluded from the ATCA list as discussed above). 
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Table 4 
VLA Radio Core Data 



PKS Name 


Other Name 


Radio Class 


z 




mjy) a 




3caml 


") 


sL G ^(mJy) 


cr 




Beam 


(") 


R 


0023-26 




CSS 


0.322 


995.8* 


7.0 





7x0 


3 


-7.1) 


820* 


9.0 





5x0 


2 


-10.2)<0.0059 


0034-01 


3C015 


FRII 


0.073 


31.5 


0.5 





4x0 


3 


51.4) 


36.9 


0.9 





3x0 


2 


53.7) 


0.0218 


0035-02 


3C17 


(FRII) 


0.220 


320.6 


1.9 





4x0 


4 


45.7) 


249.8 


4.2 





2x0 


2 


-79.2) 


0.2870 


0038+09 


3C18 


FRII 


0.188 


55.6 


0.5 





3x0 


3 


-51.1) 


68.2 


0.9 





2x0 


2 


51.4) 


0.0387 


0105-16 


3C32 


FRII 


0.400 


< 1.1 


0.4 





5x0 


3 


4.8) 


< 2.3 


0.8 





3x0 


2 


5.9) 


<0.0015 


0117-15 


3C38 


FRII 


0.565 


< 2.4 


0.8 





5x0 


3 


6.6) 


< 2.8 


0.9 





3x0 


2 


10.3) 


0.0019 


0235-19 


OD-159 


FRII 


0.620 


< 1.7 


0.6 





6x0 


3 


8.9) 


< 2.8 


0.9 





4x0 


1 


11.1) 


<0.0016 


0347+05 




FRII 


0.339 


< 2.2 


0.7 





4x0 


3 


-55.9) 


< 2.5 


0.8 





2x0 


2 


-62) 


<0.0019 


0349-27 




FRII 


0.066 


17.5 


0.4 





8x0 


3 


15.4) 


17.4 


1.1 





6x0 


2 


15.7) 


0.0095 


0404+03 


3C105 


FRII 


0.089 


19.8 


0.5 





4x0 


3 


-25.5) 


? 







3x0 


2 


63.9) 


0.0083 


0806-10 


3C195 


FRII 


0.110 


32.7 


0.4 





4x0 


3 


-1.6) 


20.7 


0.4 





3x0 


2 


-4.0) 


0.0166 


0859-25 




FRII 


0.305 


< 1.9 


0.6 





7x0 


3 


0.6) 


< 2.5 


0.8 





5x0 


2 


2.3) 


<0.0013 


0915-11 


Hydra A 


FRI 


0.054 


211.3 


2.5 





5x0 


3 


-1.8) 


184.5 


1.7 





3x0 


2 


-2.0) 


0.0146 


0945+07 


3C227 


FRII 


0.086 


13.9 


0.3 





4x0 


3 


-32.2) 


11.7 


0.3 





3x0 


1 


-23.3) 


0.0049 


1136-13 




FRII 


0.554 


234.7 


1.4 





5x0 


3 


-0.5) 


230.4 


1.6 





3x0 


2 


-1.8) 


0.1395 


1226+02 


3C273 


C/J 


0.158 


24060.0 


102.0 





4x0 


3 


-46.7) 


19460 


185 





2x0 


2 


-59.1) 


1.1884 


1306-09 




CSS 


0.464 


505.4* 


4.7 





4x0 


3 


-2.5) 


397.9* 


3.8 





3x0 


2 


2.8) 


<0.0152 


1559+02 


3C327 


FRII 


0.104 


15.2 


0.3 





4x0 


3 


44.3) 


11.2 


0.3 





3x0 


2 


53.8) 


0.0046 


1602+01 


3C327.1 


FRII 


0.462 


63.7 


0.6 





4x0 


3 


50.9) 


65.5 


0.7 





3x0 


2 


58.3) 


0.0734 


1648+05 


Here A 


FRI 


0.154 


4.5" 
























0.0004 


1949+02 


3C403 


FRII 


0.059 


20.0 


0.0 





4x0 


4 


53.6) 


8.5 


0.4 





3x0 


2 


53.6) 


0.0060 


2135-14 




FRII 


0.200 


106.3 


1.1 





5x0 


3 


-7.7) 


44.3 


1.1 





3x0 


2 


-14.1) 


0.0577 


2135-20 


OX-258 


CSS 


0.635 


416* 


3.8 





6x0 


3 


-3.1) 


199.8* 


3.3 





4x0 


2 


-7.1) 


0.0076 


2211-17 


3C444 


FRII 


0.153 


< 1.3 


0.4 





5x0 


3 


1.1) 


< 1.8 


0.6 





5x0 


29(-7.6) 


<0.0007 


2221-02 


3C445 


FRII 


0.057 


27.3 


0.4 





4x0 


3 


46.6) 


16.1 


0.5 





3x0 


2 


14.5) 


0.0096 


2314+03 


3C459 


FRII 


0.220 


< 162. 3 C 


1.0 





4x0 


3 


49.9) < 92. 6 C 


0.9 





3x0 


2 


58.5) <0.1060 



VLBI core upper limits of compact objects in the sample from Tzioumis ct al. (2002 



0023-26 


CSS 


< 20mJy at 2.3GHz 


0252-71 


CSS 


< AQmJy at 2.3GHz 


1151-34 


CSS 


< 20mJy at 8.3GHz 


1306-09 


CSS 


< 28.5mJyat 2.3GHz 


1814-63 


CSS 


87mJy b at 2.3GHz 


2135-20 


CSS 


< 11.5mJy at 5GHz 



NOTE. — Tabled Flux density results from VLA observations. (*) total flux measurement as the resolution of the observations is 
not high enough to resolve the core in these compact objects. 

a 8.4Ghz flux taken from lGizani & Leahvl d2003fl 

b Possiblc detection at the 15cr level. 

C VLBI radio maps of PKS2314+03 presented in[Th omasson ct al. (2003) resolve the central region of this object into two compo- 
nents of which it is unclear which is the true flat spectrum core. Our images cannot resolve this central structure and therefore the 
flux measurement presented here is likely to be an upper limit of the true flat spectrum core flux. 
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In a few cases, in particular those for which the 
core was not detected, we have combined both fre- 
quencies (17 and 24 GHz) in order to improve the 
sensitivity (PKS0043-42, PKS0859-25, PKS1547- 
79). Upper limits to the flux and power of the 
undetected cores, were estimated using 3ct of the 
noise level. A few objects are known to be com- 
pact steep spectrum (PKS0252-71, PKS1151-34, 
PKS1814-63, PKS1934-63). These sources are al- 
ready part of the list of ATCA calibrators and in- 
deed the resolution of our observation is not good 
enough to pinpoint a possible core. We have used 
these objects as a check of on our measured flux 
level (they agree very well with the data from the 
calibration archive and they show no indication of 
variability). For these sources, core fluxes or up- 
per limit have been m easured from VLBI maps in 
iTzioumis etafl (|2002l) (listed at the end of Table 

n. 

Finally, in a number of objects (PKS0039-44, 
PKS0043-42, PKS0409-75, PKS0442-28, PKS0859- 
25, PKS1355-41, PKS1547-79, PKS1932-46, PKS 
2250-41) we also detected radio emission from 
other regions of the radio galaxy and, in particu- 
lar, from the hot spots. Two examples of this are 
shown in Fig. [2 These maps illustrate the fact 
that at the high radio frequency and resolution of 
the observations, the maps are dominated by the 
compact cores, hotspots and knot features along 
the jets, and the diffuse lobe emission is resolved 
out. 

4.3. VLA Observations 

We observed 27 objects with the VLA at fre- 
quencies of 14.9 and 22.4 GHz. The observations 
were carried out between 4 Aug, 2-3 Sep 2006 in B- 
array configuration. Due to servicing, 10 out of 27 
dishes were not operating during the observation 
period. For one object (PKS0442-28) the quality 
of the VLA data turned out to be quite low due 
to the declination of the source, we therefore have 
excluded it from our VLA list and we will use only 
the ATCA data. 

Each source was observed for about 15 min 
per frequency, separated in three scans of 5 min 
each. In addition, a reference pointing observa- 
tion was carried out before every source. The ob- 
servations at the two frequencies were interleaved, 
to provide a better uv-coverage. As flux calibra- 
tors we used PKS0713+438 for the U-band (14.4- 



14.9GHz), and PKS1331+305 (3C286) for the K- 
band (22.1-26.0GHz). These calibrators were ob- 
served a few times during the run. Data splitting 
and opacity corrections were done inside AIPS, 
while the rest of the data reduction was done us- 
ing the MIPJAD software package. Images of the 
separate frequencies were obtained for each source 
and combined images were also made in order to 
attempt to detect those objects with faint cores, 
to no avail. Table Q] summarises the results of the 
VLA observations. 

Overall, from the VLA observations, we de- 
tected cores in 58% of the sources (15 cores de- 
tected in the 26 successfully observed). 

Again, in 7 of the VLA objects (PKS0117- 
15, PKS0347-05, PKS0859-25, PKS0915-11, PKS 
1136+07, PKS1602-01, PKS2314-03) we also de- 
tected radio emission from extended regions of the 
radio galaxy and, in particular, from the hot spots. 

4.4. Notes on the radio observations 

We chose a high resolution array configuration 
in order to try and detect the high frequency ra- 
dio cores. Due to this the resolution of the images 
resolves out the extended structure of the sources. 
Therefore it would be impossible to provide accu- 
rate measurements for the total emission from our 
sources with these observations at 15 to 22 GHz. 

Also, in the final column of the ATCA and VLA 
data tables we present estimates of the orienta- 
tion sensitive R parameter defined as S core / {Stot — 
S core ). In order to determine this we have used es- 
timates of the to tal radio flux m e asure d at 5GHz 
Stot taken from iMorganti et al.l (|19931 ). and the 



core measurement is taken from the lower of the 
two observed frequency core detections of either 
the ATCA (18 GHz) or VLA (15 GHz) observa- 
tions. Here we assume that the flat spectrum core 
flux remained constant between 5 GHz and the 
higher frequencies. The analysis of these data will 
be addressed in the detailed discussion paper to 
follow. 

5. Non-thermal Contamination 

MFIR emission can potentially have a thermal 
and/or non-thermal origin in radio loud AGN: 
thermal emission from re-radiation of AGN light 
by dust, and non-thermal synchrotron radiation. 
Therefore, in order to investigate the MFIR dust 
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Table 5 
SED Truth Table 
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Note. — Tabl({5] Truth table investigating the potential non-thermal contamination of the thermal MFIR. emission, see 
EH and 3B~3l 



emission we need to be aware of the potential con- 
tamination by non-thermal sources. There are two 
important components to consider when investi- 
gating the non-thermal contamination in a sample 
of steep-spectrum radio-loud AGN. 

• Total steep spectrum emission: Steep 
spectrum components such as lobes and 
hotspots could be potential non-thermal 
contaminants of the MFIR emission. De- 
spite the fact that these components fall in 
flux towards higher frequencies, they could 
still potentially dominate the flux of an ob- 
ject at MFIR wavelengths. However, the 
proportion of this non-thermal contamina- 
tion also depends on how much of the steep 
spectrum emission regions lies within the 
MFIR instrument beam (in this case the 
Spitzer beam). Thus, the degree of con- 
tamination by steep spectrum components 
for most of the nearby extended sources 
in our sample will be insignificant because 
they have a high proportion of their syn- 
chrotron emitting lobes/hotspots far outside 



the beam. The contamination is likely to be 
most important for compact objects such as 
CSS/GPS, and for extended FRII sources at 
high redshift that appear small because of 
their large distances. At the average redshift 
(z = 0.244) of our complete sample, objects 
with diameter < 80kpc would fit entirely 
within the 18" Spitzer beam at 70/im. 

• Flat spectrum core/jet components: 

Flat spectrum core/jet component can also 
be a potential non-thermal contaminant of 
the MFIR emission. We can detect strong 
non-thermal jet emission at optical wave- 
lengths in some objects (e.g. 3C273), so it 
is reasonable to assume that this emission 
may contribute to the observed flux at MFIR 
wavelengths. Quasars and BLRG are known 
to have str ong non-thermal bea med core 
components (jMorganti et al ] |l997f) . with the 
strength of this contamination is dependent 
on the orientation of the object to the line 
of sight. 
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Identifying the non-thermal contribution to the 
MFIR is not a trivial task, because the spectral 
shape of radio-loud AGN is often poorly sampled 
in key areas such as the sub- mm and far-infrared. 
Although some flat spectrum quasars and BL 
Lac/Blazar objects are clearly dominated by their 
non-thermal component at sub-mm and infrared 
wavelength, our steep spectrum selected sample 
has been chosen to avoid such objects, in order 
to investigate the thermal emission from dust. 
Despite this, the contribution from non-thermal 
emission from objects with relatively strong radio 
cores and from steep spectrum lobe/hotspot com- 
ponents in the Spitzer beam, remains a possibility, 
and should be carefully investigated. 

5.1. Spectral Energy Distributions 

In Figures[3]to[8]we present spectral energy dis- 
tributions for the entire sample described in sec- 
tion^ The data plotted includes the results of the 
Spitzer, ATCA and VLA observations presented in 
this paper, along with IRAS and ISO photometry 
total flux measurements and other total flux data 
taken from NASA/IPAC Extragalactic Database 
(NED). Additional core radio dat a is taken from 
Moreanti et ai] (ll993L Il997l Il999h . 

Total flux observations from ATCA of a large 
proportion o f the sample have been made at 18.5 
and 22 Ghz ( Ricci et al. 20061 ). which would com- 
plement our core measurements well. Unfortu- 
nately the atmospheric phase stability during the 
latter observations was poor, which may have con- 
tributed to the fact that, when plotted in the 
SEDs, many of the values seemed unphysical com- 
pared to the other data sets. Thus these results 
are not included in our SED plots. 

The solid line in each plot represents a single 
power-law fit to the total radio flux data of the 
objects between 10 9 to 10 10 Hz. A few objects 
have an additional dashed line fitted to the data 
between 10 9 and 10 11 Hz, when good data above 
10 10 Hz is available. Unfortunately this is a rela- 
tively under-sampled region in our data set, and 
we can fit this extra line to only 6 of the objects 
in our sample. In the following analysis this sec- 
ond fitted line, which takes better account of any 
high frequency steepening/flattening of the SED, 
is used in preference to the first solid line fit, where 
applicable. 



5.2. Lobe/hotspot synchrotron contami- 
nation 

By using the fitted power laws to extrapolate 
the total synchrotron radio emission through to 
the MFIR, it is possible to gain an indication of 
whether a contribution to the MFIR flux is pos- 
sible from the non-beamed synchrotron emitting 
lobes and hotspots. However, Spitzer 1 s beam size 
of 6, 18 and 40" at 24, 70 and 160/xm respectively, 
means that for many targets in our sample most 
of the extended steep spectrum emission lies well 
outside the beam. In Table[5]we address this issue, 
presenting a truth table from a visual analysis of 
the SEDs, where we have used the 70/mi flux as a 
reference point for the analysis. 

For an object to be considered as a candidate 
for possible steep spectrum synchrotron contami- 
nation, a source was required to fulfil two criteria. 
Firstly, the power law extrapolation of the total 
(steep spectrum dominated) radio lobe/hot-spot 
flux was required to fall close to or above the 70/im 
flux. Secondly, a substantial fraction of the total 
radio emission at 5 GHz (> 30%, as determined 
from visual inspection of radio maps) was required 
to fall within the Spitzer beam at 70/im. A tick 
in columns 5 and 11 in Table [5] indicates that an 
object fulfils both of these criteria. We find that 
only 7 out of our complete sample of 46 fulfil both 
these criteria, therefore indicating that these ob- 
jects have the potential for contamination by their 
extended synchrotron components that fall within 
the Spitzer beam. We emphasise that this is a con- 
servative estimate, given that the steep spectrum 
synchrotron component is likely to fall quicker 
than the fitted power-law at higher frequencies due 
to spectral aging of the electron population. Ev- 
idence for such spectral steepening can be seen 
in the SEDs of PKS0105-16, PKS0625-35 and 
1549—79 (S ee Figures 5 and 6). Indeed, other au- 
thors (e.g. IShi et all liobH IClearv et alJl2007h fit 
a parabola to the steep spectrum components in 
order to take this high frequency steepening into 
consideration. Therefore, our power-law extrapo- 
lation is likely to provide an over-estimate of the 
rate of steep spectrum contamination at MFIR 
wavelengths in our sample. We believe that more 
detailed parabolic or multiple power-law fits to the 
total radio flux spectra are unwarranted given the 
lack of sub-mm data for most of the objects. 
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5.3. Core/jet synchrotron contamination 

The thermal dust contribution to the MFIR ob- 
served flux can be identified by extrapolating the 
slope of the total radio flux data to higher frequen- 
cies and interpreting ever ything above that slope 
as the thermal bump (e.g. iHughes et al. 1997 ; see 
previous section IST2")) . However, this does not take 
into account any contribution one might have from 
a flatter spectrum core/jet component. 

Using our high frequency radio core observa- 
tions, and additional data from the literature, we 
have plotted the radio core fluxes and upper limits 
on the SEDs of most objects, in order to investi- 
gate their possible contribution to the MFIR flux. 
Again, using the 70/im flux point as a reference 
and assuming the core spectral shape to be flat, we 
then deem non-thermal contamination to be pos- 
sible for those objects whose high frequency radio 
core/jet fluxes lie level with or above the 70/im 
flux. 

Note again that this criterion is conservative in 
the sense that the SEDs of flat spectrum radio 
core components may not remain fiat up to MFIR 
wavelengths as we have assumed. Indeed, the 3 
flat spectrum, core dominated objects in our sam- 
ple (3C273, PKS0521-36 and PKS1549-79) all 
show significant declines between the radio and 
the sub-mm/MFIR. In the case of 3C273 this de- 
cline is 2 orders of magnitude between the radio 
and MFIR. 

We have also investigated the alternative of fit- 
ting a power-law to the core radio data and extrap- 
olating this through to infrared frequencies. Un- 
fortunately, due to the small number of core data 
points and relatively narrow wavelength range, 
there is relatively large uncertainty in the extrap- 
olations of the core data to the MFIR. Future ob- 
servations in the sub-mm region for the complete 
sample will allow us to constrain the possible con- 
tribution of the flat spectrum core-jet components 
to the MFIR more accurately. 

In columns 6 and 12 of Table [5] we present the 
results, finding that 11 out of 46 objects in our 
complete sample have a possibility of contamina- 
tion of their MFIR emission from flat spectrum 
core/jet components. 



5.4. Level of potential non-thermal con- 
tamination 

Overall, we have found that 15% of our com- 
plete sample have potential contamination from 
the steep spectrum components (lobes/hotspots), 
and 24% have potential contamination from a 
core/jet component. Out of the 46 objects in our 
full complete sample, a maximum of 30% have 
the possibility of significant non-thermal contam- 
ination of their thermal MFIR emission from ei- 
ther steep spectrum components and/or a core/jet 
component. However, as discussed above ( §5.2l fe 
§5.3|) these are likely to be conservative estimates, 
because the strength of both the flat core and 
steep spectrum components conceivably declines 
towards MFIR wavelengths faster than our simple 
analysis assumes. 

We stress to the reader that we have deliber- 
ately chosen a conservative approach to estimat- 
ing the non-thermal contamination in this study. 
The results we present here are therefore likely to 
represent an upper limit to the true degree of non- 
thermal contamination in our sample. 

5.5. Variability 

The two core-dominated objects 3C273 and 
PKS0521— 36, observed by Spitzer, were included 
in our observations and this study for comparative 
purposes, because of the large amount of previ- 
ous data available. It is possible to identify their 
non-thermal core /jet component contribution at 
MFIR wavelengths if a relatively short timescale 
flux variation is observed, since the synchrotron 
beamed component is thought to originate from a 
very compact region. 

The MFIR emission from 3C273 has been noted 
to be variable, decreasing between IRAS (1983) 
and ISO (1995-98) observations in a way that is 
consistent with monitorin g at other wavelengths 
( Meisenheimer et al.l 12001 ). We can report that 
3C273 has declined further still in MFIR emission 
since these previous observations, by a similar fac- 
tor x 2 between the IRAS and ISO data (see Figure 
O. Because of the time scale of this variation, it 
is most likely due to a change in core synchrotron 
emission, since the cool extended dust is unlikely 
to vary on the current observing timescales of 
decades. This implies that, claimed ISO detec- 
tions of thermal emission from dust underlying the 
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powerful non-thermal emission in 3C273, are un- 
likely to have been correct. The far-infrared emis- 
sion now lies well below that of the supposed ther- 
mal bump. Our data provide no evidence for the 
detection of the thermal component in the MFIR 
SED of 3C273. 

The MFIR emission of the BL Lac object 
PKS0521-36 is also known to be dominated by 
its non-thermal component, however we do not 
see such strong evidence for variability. 

5.6. Comparison with previous studies 

Previous MFIR investigations have considered 
the relative contributions of thermal and non- 
thermal emission in radio galaxies and quasars 



(Hcs ct al. 1995, van Bemmcl et al. 1998, Pollctta ct al 
2000llClearv et al.ll2007l IShi et al.ll2005h. Studies 



such as thos e bvlPolletta et al. ( 2000h . lClearv et al 
(120071) and IShi et all (|2005h fit radio to MFIR 
SEDs using various synchrotron and thermal in- 
frared emission components. Our findings agree 
with the main conclusions from these studies that 
the proportion of objects with non-thermal con- 
tamination of the MFIR by synchrotron emitting 
components is small for an unbiased sample of 
radio-loud AGN. 

In the context of unified schemes, we can also 
investigate the optical classifications of those ob- 
jects we believe have a possibility for core non- 
thermal contamination of the MFIR thermal emis- 
sion. Out of the 11 objects with possible non- 
thermal contamination from the cores, 6 are clas- 
sified as WLRG and 5 as BLRG/Quasars. A large 
proportion of BLRG/Quasars with non-thermal 
contamination is what we might expect if these 
objects have a beamed component orientated close 
to the line of sight. Overall our results our consis- 
tent with unified schemes that require a beamed 
component f or the BLRG/Quasar objects (e.g. 
lBarthel|[l989h . 

For the WLRG - some of which are FRI galax- 
ies - significant non-thermal contamination is con- 
sistent with the relatively stronger cores detected 
in FRI sources in general, as well as the detection 
of non-therm al cores in such sources at optical/IR 



wave lengths (jCapetti et al.ll2007l iChiaberge et al 
1999b . 



6. Summary 

We have presented new MFIR and high fre- 
quency radio core data for a complete sample of 
powerful southern radio galaxies. We detect ob- 
jects down to 0.5mJy at 24//m, 8.4mJy at 70/im 
and 83.6mJy at 160/xm, obtaining a uniquely high 
detection rate at 24 and 70/im. 

In addition, we have presented new high fre- 
quency radio core flux measurements, detecting 
radio cores for 59% of our complete sample. With 
these data we have made a conservative estimate 
of the non-thermal contribution to the MFIR con- 
tinuum. Careful analysis of the SEDs for our en- 
tire sample reveals that non-thermal contamina- 
tion of the MFIR is possible for a maximum of 
30% of the sources in our sample. 

An in-depth analysis of these data will be pre- 
sented in a second paper (Dicken et al. 2008, in 
preparation) . 

This work is based (in part) on observations 
made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a con- 
tract with NASA. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under 
contract with the National Aeronautics and Space 
Administration. 

D.D., acknowledges support from the STFC. 
Facilities: Spitzer {MIPS), ATCA, VLA. 
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Fig. 3. — SEDs of all objects in our sample are presented here. Upper limits are represented by downwards pointing arrows. The blue line 
is fitted to the data from 10 9 to 10 10 Hz and the red dashed line is fitted to the data points from 10 9 to 10 11 Hz. The V LA data in CSS 
object PKS0023— 26 is the total flux measurement and an upper limit on the core, plotted as an open circle, is taken from iTzioumis et al.l 
(l2Q02h . 
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Fig. 4.— SEDs cont. The ATCA data in th e SEDS of CSS object PKS0252-71 is the total flux measurement and an upper limit on the 
core, plotted as an open circle, is taken from iTzioumis et al. ( 20021 ). 
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Fig. 5.— SEDs cont. 
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Fig. 6. — SEDs cont. PKS1226+02 (3C273) clearly demonstrates its flat spectrum dominant characteristics from the fitted power law. Note 
the variability in the MF1R flux density of 3C273, where there is approximately a decade between observations of IRAS, ISO and Spitzer. 
This object is clearly in low synchrotron emission phase, see section [5751 The ATCA/VLA data in the SEDs of CSS objec ts PKS1151— 34 
andPKSf 306— 09 is the total flux measurement and an upper limits on the cores, plotted as open circles, are taken from iTzioumis et al 
(l2002h . 
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Fig. 7.— SEDs cont. ATCA/VLA core measurement in PKS1648+05 from iGizani k Leahvl (|2003f). PKS1938-15 was not observed in 
our high radio frequency program. The ATCA data in the SED of CSS obje ct PKS1814— 63 i s the total flux measurement and a possible 
detection of the core at a 15 a level, plotted as an open circle, is taken from iTzioumis et al.l ( 2002 ). We have included this PKS1938— 15 
as a candidate for steep spectrum non-thermal contamination, despite the extrapolation falling below the 70/xm. We believe that a jet 
component may be contributing to the turn up in the high frequency radio region its SED. 
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Fig. 8. — SEDs cont. The VLA data in the SEP of CSS objec t PKS2135— 20 is the total flux measurement and an upper limit on the core, 
plotted as an open circle, is taken from iTzioumis et al. ( 2002 ). 



